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Abstract 

A basic aspect of the kinetic descriptions of incompressible fluids based on an inverse 
kinetic approaches is the possibility of satisfying an H-theorem. This property is in 
fact related to the identification of the kinetic distribution function with a proba- 
bility density in an appropriate phase space. Goal of this investigation is to analyze 
the conditions of validity of an H-theorem for the inverse kinetic theory recently 
proposed by Ellero and Tessarotto [2004, 2005]. It is found that the time-dependent 
contribution to the kinetic pressure, characteristic of such a kinetic model, can al- 
ways be uniquely defined in such a way to warrant the constancy of the entropy. 

PACS: 47.27.Ak, 47.27.Eq, 47.27.Jv 
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1 Introduction 



The possibility of defining an inverse kinetic theory for the 3D incompressible 
Navier-Stokes equations (INSE; see Appendix A), recently pointed out [l[ 
^ 1^ raises the interesting question whether the relevant kinetic distribution 
function satisfies an H-theorem, namely the related (Shannon) kinetic entropy 
can be specified in such a way to result monotonically non-decreasing in time 
for arbitrary fluid fields defined in a internal domain, to be identified with 
a bounded three-dimensional domain i7 C M'^ (fluid domain). The validity 
of such a theorem is in fact a sufficient condition of strict positivity of the 
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distribution function. Therefore, this result is important in order to estabhsh 
also the consequent interpretation of the kinetic distribution function, to be 
suitably normalized, in terms of a probability density on an appropriate phase 
space r. 

Purpose of this Note is to evaluate the kinetic entropy and the related entropy 
production rate for strictly positive, suitably smooth, but otherwise arbitrary, 
distribution functions /(x,t) which correspond to an arbitrary strong solution 
of the initial-boundary value problem of INSE (see Appendix A and Refs. 0, 
1^), defined for an internal domain of M^, and for kinetic distributions functions 
which are not necessarily Maxwellian. We intend to prove that, under a suit- 
able assumptions, which involve the specification of the only non-observable 
free parameter of the theory, a strictly positive time-dependent additive con- 
tribution to the kinetic pressure, the kinetic entropy results identically con- 
served, thus yielding an H-theorem for the kinetic distribution function. The 
conclusion is obtained invoking mild assumptions both on the initial kinetic 
distribution function and on the fluid fields {p(r,t) = po > 0, V(r, t), j9(r, t)}. 
In the sequel we intend to show that the proof of this statement, and hence 
of the positivity of the kinetic distribution function, relies essentially only on 
the following assumptions: 

a) the strict positivity of the initial distribution function /(x,to) in the whole 
phase space F; 

b) the requirement that /(x,to) belongs to the functional class 



and results suitably summable in T and smooth in F x /; 

c) no-slip, Dirichlet boundary conditions are imposed on the fluid fields (see 
Appendix A); 

d) the total mass of the fluid is conserved, i.e.. 



e) the fluid is subject to a volume force density f (r,t) which is assumed suitably 
smooth, precisely at least 



/c 



(1) 




n 



f(r,t)eC(^'°)(fi X /), 
f(r, V, t) G C(°)(n X /), 

where C(^J)(fi x /) = C'^''\n) x C^^)/), with ij G N; 



(3) 
(4) 
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f ) the existence of a strong solution of INSE exists in Q x 7 which belongs to 
the functional class: 



< 



V{r,t),p{r,t) E C(°)(^ X I), 
V{r,t),p{r,t) e C(2,i)(f] X /), 



(5) 



with I denoting the time axis, generally to be identified with a finite subset 



g) the determination of the time-dependent term of the kinetic pressure (pi). 
It if found that pi can always be defined in the same time interval 7, up to 
an arbitrary positive constant, in such a way that the entropy production rate 
vanishes identically. 

The basic result can be summarized as follows. 

Theorem 1 - H-theorem (entropy conservation for the inverse ki- 
netic theory of INSE) 

Let us denote by f{'x.{t),t) the kinetic distribution function solution, assumed 
to exist and result suitably regular, of the inverse kinetic equation represented 
in the form 



X = (r, v) being a state vector belonging to the phase space T = fl x U and 
U — M.^ the velocity space. /(xo,to), the initial kinetic distribution function, 
is assumed to be suitably smooth in the sense of Eq. (1), strictly positive and 

summahle in T for appropriate weight functions. Moreover, Tt^to suitably- 
defined diffeomorphism, denoted as Navier-Stokes evolution operator, such that 
Vt G /, including the initial time to G /, the fluid fields at any time t E I, 
{po,'V{r,t),p{r,t)} are uniquely defined by the following velocity moments of 
the kinetic distribution function: 



of R; 




(6) 



Po 



j d^vf{x,t) = Po 



I d^vTt,tJ{^o,to), 



(7) 



V(r,t) 



1 



/ 



d^vvf{x,t) 



1 



/ 






where 



pi{r,t) = I dvEf{x,t) = I dvETt,tJ{^o:to). (10) 



Here E = ^u"^, with u = v — V, while Pg = Po{t) is a smooth strictly positive 
function. Then, the following results follow: 
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1) provided the initial distribution function /(xo,to) results strictly positive 
and suitably smooth in the sense of (1), it follows that the velocity moments 

Fciv^t) = / (i^t>G(x,t)/(x, t), where G{yi,t) = l,v,E = ^u'^, uu,uE and 

u _ 
u = V — V, exist, are continuous Q x I and suitably smooth in Q x I in the 
sense of the settings (2), (3), (4) and (5). In addition, the kinetic entropy 

S{t) = - y'cix/(x,t)ln/(x,t) (11) 
r 

exists and is suitably smooth in I ; 

2) provided the kinetic pressure p{r,t) = J d^vEf{:K,t) is suitably prescribed, 

u 

there results identically in I 

S{t) = S{to) (12) 

(law of entropy conservation). 

For greater clarity, in Sec. 2 the relevant aspects of the inverse kinetic theory 
previously developed are recalled j^. This is useful to introduce the Navier- 
Stoker evolution operator T^^to ci-nd define the related probability density in the 
phase space T. As a consequence, it is immediate to prove that Tf^to conserves 
probability in the same space. Subsequently, in Sec. 3, the kinetic entropy 
the S{t) and its time derivative dS{t)/dt are evaluated. It is found that by 
suitably defining the non- negative kinetic pressure (pi) the entropy production 
rate can always be set equal to zero in the whole time interval I in which by 
assumption a strong solution of INSE exists. Implications of the result are 
pointed out. 



2 The Navier-Stokes evolution operator 



In this Section we briefly recall the formulation of the inverse kinetic theory 
developed in . This is useful in order to identify the Navier-Stokes evolution 
operator Tf which determines the evolution of the relevant kinetic distribu- 
tion function /(x, t) in the extended phase space F x / ( where T = fl x U , 
with f/ = denoting a suitable "velocity" space). 

The result is obtained by requiring the /(x, t) obeys a Vlasov kinetic equation 
of the form 

L(F)/ = 0, (13) 
where x = (r, v) is a state vector spanning the phase space F, L the streaming 
operator 
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and X a vector field of tfie form X(x,t)= {v, F(x, t)} . This equation can for- 
mally be written in integral form by introducing the initial kinetic distribution 
function /(xo,to) = /^(xo), defined at the initial time to G /, and the flow, 
i.e., the the diffeomorphism Xo x(t) =Tt^to^o generated by the vector field 
X, via the initial-value problem 

|x = X(x,t) 

dt V > ; ^^^^ 

x(to) = Xo 

and its related evolution operator Tt^to {Navier- Stokes evolution operator). 
This implies that provided the solution of the initial-value problem exists, is 
unique and suitably smooth, the Jacobian of the flow Xq x(t), J(x(t),t) = 
is non singular and reads 

J(x(t),t) =exp|^*dt'^-F(x(0,0|. (16) 

Thus, the evolution operator, acting on the kinetic distribution function /o(x^) 
results 



/(x(t),t) = r,,Jo(xJ = /o(xJexp|-£^;t'^.F(x(t')/ 



(17) 



Previously it has been shown [4] that the functional form of the vector field 
X and of the "mean- field force" F(x, t) yielding and inverse kinetic theory for 
INSE can be uniquely determined under suitable prescriptions. These include, 
in particular, the requirements that: 

1) the local Maxwellian kinetic distribution function 

Po 



/M(x,t;V,pi) = ^— expj-X^I 



th 



[where X"^ = v^f^ = 2pi/ po and u is the relative velocity u=v — V(r,t)] re- 
sults a particular solution of the inverse kinetic equation if and only if {p, V,p} 
satisfy INSE; 

2) suitable bounce-back boundary conditions are imposed for the kinetic dis- 
tribution function on the boundary 6Q 



3) the moment equations corresponding to the velocity moments MG{r,t) = 
J (PvG{'x, t)/(x,t) for G(x, t) = 1, v/po, E = coincide with the differential 
equations of INSE; 
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4) the fluid fields po and V(r, t) are identified respectively with the velocity 
moments for G{^,t) = 1, v/p^; similarly, the fluid pressure p{r,t) is defined in 
terms of the kinetic pressure Pi(r,t) [see Eq.(9)] by requiring 

p{r,t)^p,{r,t)-Po. (19) 

It is obvious, in order that Vp = Vpi, that Pg can be in principle an arbitrary 
strictly positive function independent of r. Thus it can always be assumed to 
be Vt e / a smooth function of t. The resulting form of F(r, v,t) is recalled in 
Appendix B. It implies 



2pi 

1 



2p 

where = ^ + V-^ is the Lagrangian (or convective) derivative. For f = /m 
it becomes 



V n 



qI + — u- vn 

J pi — 



^ ■ F(x,t) + -u-Vp, (21) 

c/v 2pi Dt pi 

which implies that if / = /m at a given time These expressions, in particular 

(20) and (21), permit to determine uniquely the Jacobian J(x(t),t) and the 

evolution operator T^^to- Thus, introducing the normalized kinetic distribution 

function 

/(x,t) = l/(x,t) (22) 

Po 

and requiring that the initial kinetic distribution function f{x.o,to) = foi^o) 
results at least of class C*^^^ (F x /) and summable in F it follows 

dx(i)/(x(0, t) = d^{to)fiyiito),to) = d^ofoM, (23) 
and in particular 

J dK{t)f{K{t),t) = J dxj„(x„) = 1. (24) 



In order to prove that f{x{t),t) can be interpreted as probability density in 
the next section we intend to establish an H-theorem 



3 Shannon kinetic entropy 



In terms of the Navier-Stokes evolution operator and Eq.(20) [or (21) in 
the case in which /(x,t) coincides with a local Maxwellian distribution (18)] 
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it is now immediate to evaluate the Shannon kinetic entropy associated to the 
kinetic distribution function /(x.t), namely S{t) = — /p (ix/(x,t) In /(x,t). 
Let us assume, for this purpose that the initial kinetic distribution function 
/(xq, to) = /o(xo) be defined in such a way that it results strictly positive in 
r, at least of class C^^^(r x I) and summable in F so that the Shannon kinetic 
entropy S{to) = — /p (ixo/(xo, to) In /(xo, to) results defined and at least of 
class C^-^\l). Thanks to the integral kinetic equation (17) and the condition 
of conservation (23) it follows that S{t) and S{to) are elated by means of the 
equation: 

S{t) = S{to) + d^oU^o) [ dt'-^ ■ F(x(t'),i')- (25) 

Therefore, the entropy production rate -§iS{t) results 

|5(t) = /^dx/(x,t)|--F(x,t), (26) 

where ^ ■ F(x,t) is given or an arbitrary kinetic distribution function by 
Eq.(20) [or (21) for the Maxwellian case]. It follows 



d_ 
di 



V n 



Q 



(27) 
(28) 



Since Po + p{r,t) and /(x, t) are strictly positive, we can always define Po{t) 
so that in the finite time interval / there results identically 



-Indr 



n P„(t)+p(r,t) 



V n 



Q 



(29) 



■/n^^Po(t)+p(r,t) 

In the case in which the initial condition /o(xo) coincides identically with /m 
in the whole phase space F it follows in particular 

-Indr ' ^ 



d_ 
di 



Po(t)+p{r,t) Dt 



P 



Indr 



1 



(30) 



Po(t)+p{r,t) 



Therefore, condition (29) [or (30) in the Maxwellian case] implies that in the 
same time interval / the entropy production rate must vanish identically, i.e.. 



d_ 

di 



S{t) = 0. 



(31) 



We stress that Eq.(31) holds, in principle, for an arbitrary initial condition 
Po{to) = Poo > with Poo suitably large. Therefore, the kinetic pressure pi. 
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given by Eqs.(9) and (10), remains still non-unique since it is determined in 
terms of Eq.(29) only up to an arbitrary positive constant Poo- 
It follows that the Shannon entropy for the kinetic distribution function /(x,t) 
results always conserved by imposing a suitable prescription on the the time- 
dependent part of the kinetic pressure Po(t). Since the latter is unrelated to 
the physical observables (i.e., the fluid fields) the constraint condition imposed 
on the kinetic pressure [respectively (29) or (30)] can always be satisfied. As 
a consequence, with such prescriptions the normalized kinetic distribution 
function /(x, t) can be interpreted as probability density. 



3. 1 Conclusions 



In this paper the condition of positivity of the kinetic distribution function 
f{^,t) which characterizes the inverse kinetic theory recently developed for 
the incompressible Navier-Stokes equations has been investigated j3, 0|- We 
have proven that the Shannon entropy is exactly conserved for arbitrary ki- 
netic distribution function, provided the kinetic pressure is suitably defined 
and the initial kinetic distribution function results positive definite and suit- 
ably regular. As indicated, these conditions can always be satisfied without 
imposing any constraint on the physical observables, here represented by the 
fluid fields {po, V,p} . 

The conclusion applies in principle to arbitrary, suitable smooth in the sense 
(5), strong solutions of INSE which are defined in three dimensional internal 
domains of M^. Assuming, mass conservation and no-slip boundary conditions 
(i.e., Dirichlet boundary conditions) on the boundary 6Q, the same result 
holds also for non-isolated systems characterized by moving boundaries. In 
addition, arbitrary volume forces which satisfy (3), (4) or analogous surface 
forces obtained by applying a non-uniform pressure on the boundary 6Q, can 
be included. 

An immediate consequence of the H-theorem here obtained is the possibility 
of imposing the maximum entropy principle in order to determine the ini- 
tial kinetic distribution function /o(x), i.e., requiring the variational equation 
SS{fo) = subject to suitable constraint equations. Thus, for example, the lo- 
cal Maxwellian distribution (18) is obtained by imposing solely the constraints 
provided by the moments (7), (8), (9) and (10), to be considered as prescribed. 
However, in principle, the variationalprinciple can also be used to determine 
non-Maxwellian initial distributions jsf. 



These results appear significant both from the mathematical viewpoint and 
the physical interpretation of the theory. 
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4 Appendix A: INSE 



The incompressible Navier-Stokes equations (INSE) are defined by the follow- 
ing set of PDE's and inequalities for the fluid fields {p, V,p} 



-p + V-(pV) = 0, (32) 

p-^ V + Vp + f - /iV V = 0, (33) 

V-V = 0, (34) 

p(r,t)>0, (35) 

p(r,t)>0, (36) 

p(r,t)=p, >0. (37) 



The first three equations (32), (33) and (34), denoting respectively the continu- 
ity, forced Navier-Stokes and isochoricity equations, are assumed to be satisfied 
in the open three-dimensional set f2 C (fluid domain) and in a possibly 
bounded time interval / C R, while the last three inequalities, (35)-(37) ap- 
ply also in the closure of the fluid domain = U 6Q. Here the notation is 
standard^. Hence = ^ + ^ ■ V and p = vpo > is the constant fluid vis- 
cosity, with u the related kinematic viscosity. The volume force density f(r,t) 
acting on the fluid element by assumption is taken in the functional setting 
(3), (4) and (5). Consequently, the fluid flelds {V(r,t),p(r,t)} are required to 
satisfy the regularity conditions (5). The initial-boundary value problem for 
INSE is deflned as follows. The initial condition is deflned by imposing 



p{r,to) = po>0 (38) 

p{r,to)=Poir), (39) 

V(r,t,) = V,(r), (40) 

where {Vo(r),po(r)} belong to the functional class 



V„(r),p,(r)GC(o)(l]), 
V,(r),p,(r)GC(2)(fi), 



(41) 
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and moreover satisfy respectively the isochoricity condition (34) and the Pois- 
son equation 

VVo(r) = -V ■ {poV„-VV„ + f (r^)} . (42) 
The boundary conditions can be specified, for example, by means of the Dirich- 
let boundary conditions (which for the velocity are usually denoted as no-slip 
boundary conditions), i.e., letting Vt G / and imposing in each point r^y of 
the the boundary SQ 

p(-,t) = p„>0, (43) 
p{;t)=pw{;t), (44) 
V{;t) = Vw{;t). (45) 

Here {Vvi/(r,i),Pvy(r,t)} denote respectively the velocity and the pressure at 
an arbitrary point rw belonging to the boundary 6Q. fields, both required to 
belong to the same functional class (5). 



5 Appendix B: mean- field force 



For a generic (i.e., non-Maxwellian) distribution function /(x,t), the mean- 
field force F reads F(x,t) = Fo(x,t) + Fi(x,t), where Fq and Fi are the vector 
fields: 



Fo(x,t) = — [V-n - Vj9i - fl + a+i/V^V, (46) 



+ 



Po 

F 

1 



V n 




2pi L =J --J 2pi = I 2 
where = and vff^ = 2pi/po. Here Q and n are the velocity- moments 
Q = / d^vu^f, n = / d^vuuf, while f denotes the volume force density 
acting on the fiuid element and finally z/ > is the constant kinematic vis- 
cosity. In particular, for the Maxwellian kinetic equilibrium (18) there results 
n = pil, Q = 0. Moreover, a is the convective term which according to Ref. 
is uniquely defined and reads a =|u ■ VV + |VV • u. 
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